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ABSTRACT 

We compare the evolution of spectral shape with luminosity in GRS 1915+105 
with that of 'normal' black holes. The pathological variability of GRS 1915+105, which 
probably indicates a disc instability, does not require that GRS 1915+105 belongs in a 
different class to all the other objects. At comparable fractions of Eddington luminosity 
its spectra and (more importantly) apparent disc stability, are both similar to that 
seen in the 'normal' black holes. Its unique limit-cycle variability only appears when 
it radiates at uniquely high (super-Eddington) luminosities. 
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1 INTRODUCTION 

GRS 1915+105 is a spectacularly variable accreting black 
hole in our Galaxy. It was first detected in X-rays by the 
WATCH all sky monitor on board the GRANAT satellite 
in 1992 (Castro- Tirado, Brandt & Lund 1992), but achieved 
fame when radio observations showed it to be the first galac- 
tic superluminal source, implying relativistic plasma ejec- 
tion (Mirabel & Rodriguez 1994). Several more superlumi- 
nal jet sources have since been found (GRO J1655-40, XTE 
J1748-288, V4641 Sgr), while a much larger number of black 
holes (the microquasars) show jet-like radio morphology (e.g. 
Mirabel & Rodriguez 1999) . Mildly relativistic jets are now 
recognized as a common feature of black holes at high mass 
accretion rates (Fender & Kuulkers 2001). 

However, GRS 1915+105 is still unique in its variabihty 
properties. On long timescales, its outburst duration of over 
10 years is unprecedented, and completely at odds with the 
predictions of the disc instability models which can success- 
fully describe some aspects of the other black hole transient 
outbursts (King & Ritter 1998; Dubus, Hameury & Lasota 
2001; Esin, Lasota & Hynes 2000). On short timescales, the 
X-ray variability is even more singular, showing episodes 
where it continually switches between states in a quasi- 
regular way (Greiner, Morgan & Remillard 1996; Chen, 
Swank & Taam 1997). Belloni et al. (1997a) showed that 
these rapid changes were associated with the accretion disc 
spectrum switching from hot and bright, implying a small 
inner disc radius, to cooler and dimmer, with a larger in- 
ferred radius. They interpreted this as the result of a limit- 
cycle instability in the inner accretion disc, such that it is 
continually emptying and refilling. 

Following Belloni et al. (1997a, b), most models for the 



origin of this unique limit-cycle behaviour of the accretion 
disc have concentrated on the radiation pressure instability 
in a standard Shakura-Sunyaev accretion disc (Shakura & 
Sunyaev 1973) . In the a disc prescription the viscous heating 
is proportional to the total pressure Ptot = Pgas + frad, i.e. 
the sum of the gas and radiation pressures. Where radiation 
pressure dominates the heating rate is dramatically sensi- 
tive to temperature (Prad oc T* compared to Pgas oc T), and 
the radiative and convective cooling in the disc cannot keep 
pace. By itself, this would just lead to the classic thermal- 
viscous instability (Lightman & Eardley 1974; Shakura & 
Sunyaev 1976). However, there is another stable disc so- 
lution at high temperatures, where optically thick advec- 
tive cooling becomes important (Abramowicz et al. 1988). 
Putting all this together at a given radius and a gives rise 
to an S-curve on a plot of disc surface density, S versus 
mass accretion rate through the disc, m cx T'^. The lower, 
middle and upper branches of the S corresponding to heat- 
ing oc Pgas balanced by radiative cooling (stable), heating 
tx Prad balanced by radiative cooling (unstable) and heating 
oc Prad balanced by advective cooling (stable; Abramowicz 
et al. 1988; Chen et al. 1995). No steady state equilibrium 
solution is possible at mass accretion rates corresponding to 
the middle branch. Instead there is a continual limit-cycle 
with the material switching between the two stable states, 
one which has too high a mass accretion rate, leading to the 
emptying of the disc, and one with too low as mass accretion 
rate, so the disc refills. 

However, it is difficult to understand why the radia- 
tion pressure instability is not apparent in other black hole 
X-ray transients (see e.g. Gierlihski & Done 2003b, here- 
after GD03) . The only features distinguishing black holes are 
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mass and spin. GRS 1915+105 is the most massive galactic 
black hole known at '^14 M© (Greiner, Cuby & McCaugh- 
rean 2001a), but this is not significantly larger than the typi- 
cal WMq of the other black holes (Bailyn et al. 1998). It may 
well be spinning close to maximal, but so may GRO J1655- 
40 (Zhang, Cui & Chen 1997; Cui, Zhang & Chen 1998). 
If the difference is not in the black hole itself, then it must 
be connected to the accretion flow. GRS 1915+105 accretes 
at high fractions of the Eddington limit, but so do many 
other black hole transients (see e.g the review by Tanaka & 
Lewin 1995). GRS 1915+105 has a jet, which is clearly linked 
to hmit-cycle behaviour of the disc (Pooley & Fender 1997; 
Eikenberry et al. 1998) , but the ratio of radio to X-ray power 
is similar to that of other X-ray transients (Fender & Kuulk- 
ers 2001). Why then is the variability from GRS 1915+105 
unique? 

Recently, Done & Giorliiiski (2003; hereafter DG03) 
used the huge Rossi X-ray Timing Explorer (RXTE) 
database to systematically analyze the spectra of many black 
hole systems. These were all consistent with the same spec- 
tral evolution as a function of Eddington luminosity^ frac- 
tion, L/L-Edd, so these form a 'normal' black hole sample 
for comparison with GRS 1915+105. Here we analyze the 
RXTE data on GRS 1915+105 and plot it together with the 
'normal' black holes. We show that GRS 1915+105 is unique 
in being the only black hole binary which spends any con- 
siderable time above L/L^dd ^ 1, and speculate that this is 
the trigger for the limit-cycle variability in GRS 1915+105. 

The reason GRS 1915+105 alone gets to such consis- 
tently high mass accretion rates can be simply related to 
its evolutionary state. The secondary is a giant (Greiner et 
al. 2001b), so Roche lobe overflow occurs in a much wider 
binary than for a main sequence star. The 33.5-day period of 
GRS 1915+105 is by far the longest of any low-mass X-ray 
binary (LMXB) known. This implies a huge disc size, as the 
outer radius is set by tidal truncation. During quiescence the 
disc builds up an immense reservoir of mass, which can fall 
onto the black hole when the outburst is triggered (King & 
Ritter 1998). The underlying cause of all the unique long and 
short-term variability of GRS 1915+105 is then ultimately 
linked to the evolution of the huge disc structure, which 
can contain enough material to maintain super-Eddington 
accretion rates over timescales of tens of years. 



2 DATA SELECTION AND REDUCTION 

RXTE has been in operation since Jarmary 1996. The Pro- 
portional Counter Array (PCA) on board RXTE had un- 
dergone several major gain changes, marking five instrument 
epochs. Here we analyze all PCA data before the end of PCA 
epoch 4, i.e. observations between 1996 April 6 and 2000 
May 11 (epoch 5 is marked by the loss of the propane layer 
for detector 0). We use ftools 5.2 and follow standard data 
reduction procedure as recommended by the RXTE team. 
To maximize the number of observations which can be used 
we select data from PCA detector only, as it is switched 
on most of the time. The source is dramatically variable 

1 We use Eddington luminosity Lsdd = 1-26 X lO^^M/M© erg 
for a mass M throughout this paper. 



so we extract data in 128-s bins, giving a quite overwhelm- 
ing number of 16489 individual spectra. There is variability 
on shorter timescales (see e.g. Belloni et al. 2000; hereafter 
BOO) but this is the shortest timescale on which the spectra 
are generally limited by systematic errors (set to 1 per cent) 
rather than statistical errors below 10 keV. We also accu- 
mulate both light curves and spectra for each of the 500 
separate pointings, designated by unique observation IDs, 
typically each giving a few kiloseconds of data. The light 
curves for each observation ID are extracted over 3-20 keV 
with the 16-second time resolution of the Standard 2 data. 



3 SPECTRAL MODELLING 

We use XSPEC (Arnaud 1996) to fit the PCA data in 3- 
20 keV range with 1 per cent error added in each energy 
channel to represent the systematic uncertainties of the de- 
tector response. We assume the abundances of Anders & 
Ebihara (1982) unless stated otherwise, and fix the inter- 
stellar galactic absorbing column at 4.7 x 10^^ cm~^ (Chaty 
et al. 1996). For calculating luminosities we assume distance 
of 12.5 kpc (Mirabel & Rodriguez 1994 but see Fender et al. 
1999). All parameter ranges are quoted for Ax^ ~ 2.7. 

Following DG03 we fit a physically motivated contin- 
uum model consisting of a multicolour disc spectrum which 
also forms the seed photons for Compton scattering to higher 
energies. We use the diskbb model (Mistuda et al. 1984) to 
describe the disc emission. The temperature is constrained 
to fcTdisc > 0.4 keV. so as not to produce an unphysically 
large component outside of the energy range of the PCA 
data. For Comptonization wo use thcomp (not included in 
the standard distribution of xspec, see Zdziarski, Johnson 
& Magdziarz 1996), parameterized by asymptotic photon 
spectral index (F > 1) and electron temperature (fcTcomp < 
100 keV). 



3.1 Simple model 

The physically motivated continuum model does not fully 
describe the observed spectra from X-ray binaries as there is 
also Compton reflection from the disc. Following DG03 we 
approximate this by adding a broad Gaussian line (width 
fixed at 0.5 keV) and smeared edge (index for photoelectric 
cross-section fixed at -2.67 and smearing width of 7 keV; see 
Ebisawa 1991). This model spectrum is described as Ml in 
Table 2 and throughout the text. 

We use this to fit each 128-s spectrum from 
GRS 1915+105. It gives a good description of most of the 
data; only 5 per cent have reduced /v > 1.5 (where i/ 
is number of degrees of freedom). We follow the approach 
of DG03 in compressing all the resulting spectral informa- 
tion from the model fits into two intrinsic colours, which are 
calculated by integrating the (absorption corrected) model 
flux over 4 energy bands (3-4, 4-6.4, 6.4-9.7, 9.7-16 keV) 
to form soft and hard colours. These roughly describe the 
unabsorbed spectral slope from 3-6.4 keV and 6.4-16 keV, 
respectively. 

This method contrasts with that of the more usual in- 
strument colours, which are defined using the ratios of ob- 
served counts in given energy bands (as opposed to modelled 



© 2003 RAS, MNRAS 000, 1-12 



3 





1 1.5 2 0.1 0.2 0.3 0.4 

Soft colour HR2 

Figure 1. (a) shows the intrinsic colour-colour diagram for each 128-s spectrum of GRS 1915+105 using the data selection of BOO. 
These colours are corrected for interstellar absorption and for the instrument response, giving a measure of the intrinsic spectral slope in 
the 3—6.7 keV (soft colour) and 6.7—16.0 keV (hard colour) bandpass (see DG03). (b) shows the corresponding instrument colour-colour 
diagram as used by BOO. Regions marked in the same colours in panels (a) and (b) correspond to the same data. The shading indicates 
the spectral states A, B and C defined by BOO, while the solid/dotted line is again from BOO, indicating the path of a typical instability. 
The BOO data are on 1-s resolution, and so their colours can follow the very rapid variability in the flares (indicated by the dotted 
line) which are not sampled by the 128-s spectra used here. Panels (c) and (d) reproduce panels (a) and (b) but for all observations 
from PCA epochs 2-4. The intrinsic diagram looks very similar, but the instrument diagram is very different due to gain changes in 
the detector. There is no longer a unique correspondence between the intrinsic spectral shape and instrument count ratios (e.g. the red 
region indicating a particular spectral shape is split into three separate regions in instrument counts). 



flux ratios). Both methods have their advantages and disad- 
vantages. The advantage of our method is that modelhng 
corrects for both instrument response and absorption in a 
single step, and gives colours which reflect the underlying 
source spectrum. Thus many objects, with data taken from 
many different instruments and/or gain epochs, with very 
different absorption columns, can be easily plotted together 
on the same diagram. Instrument colours can be modified 
to correct for different gain epochs (e.g. Di Salvo et al. 2003, 
van Straaten et al. 2003), but different absorptions prevent 
a straightforward comparison of many different objects. An- 
other advantage of our approach is that the model spec- 



tra allow us to estimate the bolometric luminosity, which is 
surely a major physical driver of the spectral evolution. The 
disadvantage of course is that these luminosities are model 
dependent (though the colours are not as long as absorp- 
tion does not substantially affect the spectrum above 3 keV: 
DG03). More importantly, the requirement for good signal- 
to-noise spectra for the model fitting means that the intrinsic 
colours have lower time resolution than the instrument ones. 

Belloni et al. (2000: hereafter BOO) used instru- 
ment colours for their major study of the variability of 
GRS 1915-1-105 as seen in Epoch 3. Fig. la and b shows 
a comparison of our derived intrinsic colours (defined by 
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model flux ratios) for the Epoch 3 data with those of the in- 
strument colours (defined by observed count ratios) as used 
by BOO. Data from different regions of Fig. la are colour- 
coded, so that their corresponding position on Fig. lb can 
easily be identified, so that we can use the BOO insights into 
the nature of the variability as a function of colour. In par- 
ticular, they unified the incredibly rich variability patterns 
into transitions between 3 main spectral states - A, B and 
C which they associate with a disc instability (Belloni et 
al. 1997a, b). These states are shown by the shaded areas 
in Fig. lb, while the thick black line shows the track in 
colour-colour space of a dramatic transition (A-type vari- 
ability, taken from observation ID 10408-01-38-00 as shown 
in figure 2m of BOO) . The dotted line shows the continuation 
of this track which is seen in high time resolution data. BOO 
are able to get a good instrument colour on 1-s time resolu- 
tion, as opposed to the 128-s spectra needed for the model 
fitting used here. 

While we cannot follow the most rapid variability, 
Figs. Ic and Id show some of the power of our intrinsic 
colour approach when combining all the data from Epochs 
2-4. The intrinsic colour-colour diagram is very similar to 
that derived just from BOO data (compare panels a and c in 
Fig. 1), whereas the straightforward instrument colours are 
strongly affected, with data from different Epochs ending 
up at very different position in the diagram (Epoch 2 has 
much higher HRl and HR2 than epoch 3, while Epoch 4 
has much lower HRl and HR2). While these gain change ef- 
fects can be corrected for (di Salvo et al. 2003; van Straaten 
et al. 2003), the key issue in this paper is that the unab- 
sorbed flux ratios shown in Fig. Ic can be directly compared 
with the unabsorbed colours from many other black holes 
given in e.g. DG03. We can use this straightforward compar- 
ison of multiple objects to search for the differences between 
GRS 1915-1-105 and the 'normal' black holes which underlay 
its unique variability. 

3.2 Variability as a function of colour 

BOO show that the unique limit-cycle variability of 
GRS 1915-1-105 occurs in only ~50 per cent of the data. 
Rather than do a detailed classification of the variability 
pattern (as in BOO), we simply parameterize the variabil- 
ity in a given observation (one individual RXTE pointing 
with unique observation ID) by calculating the fractional 
r.m.s. (variance divided by mean intensity) from its light 
curve with resolution of 16 s. Some of the observations are 
too short for this to be reliable, so we restrict this to light 
curves in which there are more than 128 data points (i.e. 
over 2048 s of data). This fractional r.m.s. represents the 
integral of the power spectrum from frequencies lower than 
1/2048 to 1/32 Hz. The low frequency power spectrum is 
generally flatter than (e.g. Morgan et al. 1997) so the 
fractional r.m.s. is not strongly dependent on the length of 
the light curve. 

Assuming that the power spectrum is stable for each 
type of variability, we can use the individual observa- 
tion IDs identifled by BOO as being characteristic of each 
class to get an estimate of their fractional r.m.s.. The 
state transition classes ^, 8, A, k, p, v, a, /3, S and 7 have 
17, 36, 52, 54, 35, 47, 27, 49, 16 and 7 per cent variability, re- 
spectively, while the single state classes (j) and x have 7 and 



1.2 



1 - 



0.8 
o 

o 

1 0.6 



0.4 



1 ' \ < 1 < 1 

0.01 - 0.05 

0.05 - 0.15 1 
0.15 - 0.40 
0.40 - 0.72 ••. 

, . • 

•• . . • 

* «. •••• • 

. •>^>;.- 

• •• •: •• 

• * ••* 'jlk' • * 



1.2 1.4 1.6 1.8 2 
Soft colour 

Figure 2. The intrinsic colour-colour diagram of the individual 
PCA pointings with > 2048 s of data. These cover the full range 
of spectral behaviour seen in Fig. Ic. The points are colour coded 
according to the amount of fractional r.m.s. variability present in 
the 3—20 keV light curve for each observation as given in the key 
in the top-left corner of the figure. The numbered stars show the 
positions of the individual spectra used for spectral analysis in 
Sec. 3.3. 



3 per cent, respectively. While it is plain that our crude 
fractional r.m.s. criteria cannot make subtle distinctions be- 
tween the types of state transition, nevertheless it is clear 
that restricting the fractional r.m.s. to less than 5 per cent 
will select typically only the x stable state (through may 
include a little contamination by the (j> and 7 type of vari- 
ability), while selecting fractional r.m.s. greater than 40 per 
cent will pick out the dramatic state transition variability 
of type A, K, v, (3 (perhaps with a little contribution from 6 
and p states. It is these dramatic variability patterns which 
are unique to GRS 1915-1-105 which we want to identify. 

The intrinsic colours for the resulting 292 spectra were 
found using the model described above. Obviously, some of 
the spectra are averaged over periods of high variability, 
including the limit-cycle spectral transitions. These do not 
represent real physical spectra: the average of two black- 
bodies of different temperatures is not a blackbody at the 
average temperature, nor is the average of two power laws 
of different spectral slopes given by a single power law of 
the average slope! Instead, the colours extracted from these 
spectra give the time averaged position of the source on the 
colour-colour diagram. Fig. 2 shows these intrinsic colours, 
with r.m.s. variability indicated. This illustrates the result 
of BOO that the limit-cycle variability which is unique to 
GRS 1915-1-105 is associated with very particular spectral 
shapes which mark out the track of the instability. 



3.3 Detailed analysis of selected spectra 

We use Fig. 2 to select seven observations which do not show 
signiflcant variability and cover full range of behaviour in the 
colour-colour diagram. The chosen points are marked with 
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Model XSPEC components 



Description 



Ml 


wabs*smedge*(diskbb+thcomp+gaussian) 


bimplc model oi multicolour disc with thermal Coiiiptoiiization. He- 
flection features are approximated by broad Gaussian and smeared 
edge. Absorption with solar abundances fixed at Nu = 4.7 x 10^^ 

_9 

cm . 


M2 


WARS* (diskbb+thcomp+relrep) 


Multicolour disc, thermal Comptonization and its ionized Compton 
reflection, relativistically smeared. Absorption with solar abundances 
fixed at Nh = 4.7 x 10^2 cm'^. 


M3 


VARABS* (diskbb+thcomp+relrep) 


The same as M2, but with anomalous abundance absorption of Si and 
Fe. 


M4 


varabs*(diskbb+thcomp+relrep+6x Gaussian) 


The same as M3, but with narrow resonant absorption Lyman a, fi 
and 7 lines from He-like and H-like iron. The absorption lines are 
assumed to be narrow and the fija and 7/a ratios are fixed to the 
oscillator strength ratios. This is our best physical model. 


M5 


varabs*smedge*(diskbb+thcomp+gaussian) 


The same as Ml, but with anomalous abundance absorption of Si and 
Fe. This is a very fast model which gives very similaj: colours to M4. 


Table 2. Summajry of spectral models used in this paper. 




Obs. fcTdiac r 
(keV) 


^Tcomp ^disc .t/comp X 

(keV) (lO^Sergs-i) 



1 


< 0.74 


2 00+"-"5 


> 10 


2.6+1.6 


5 1+2-3 
^^•^-1.3 


16.4/36 


2 


0.821-11; 


2.15 ± 0.16 


> 7 


3.5 + 0.4 


4 2+1-5 


10.5/36 


3 


i-36l°;« 


1 97+018 

-0.85 


5 2+1-1 


< 6.3 




10.2/36 


4 


"•'^-0.38 




5 0+°-8 


< 4.8 


5.8l?1 


18.7/36 


5 


< 0.87 


r, oc--l-0.ll 
^■"'''-0.14 


5 5+1-1 


< 6.4 


6 4+^^-1 


10.8/36 


6 


qO+^-l^* 


^- " -0,15 


3.2 ± 0.2 


< 3.5 


11-1-1;? 


5.1/.36 




1 IV I"-"- 


< i.2() 


0,U2 


(i,()l ± 0,05 


n.')s ± 0.01 


02,2/36 



Table 3. Best-fitting parameters for the selected spectra using the standard model Ml (see Table 2). Luminosities have been calculated 
for the distance of 12.5 kpc and the disc inclination of 66°. Eddington luminosity for a 14 M© black hole is 17.6 x 10^8 erg s~i. 



Obs. 


Obs. ID 


Live time 


Count rate 


r. m. s. 






(s) 


(s-i) 


per cent 


1 


30703-01-34-00 


2400 


873+2 


4.0+0.3 


2 


20402-01-29-00 


4208 


1157+3 


2.9+0.2 


3 


10408-01-27-00 


7552 


1562+3 


1.9+0.1 


4 


30402-01-09-00 


3712 


1751+4 


1.5+0.1 


5 


30184-01-01-00 


15184 


1896+4 


1.3+0.1 


6 


20402-01-55-00 


7552 


3267+8 


3.9+0.2 


7 


30703-01-08-00 


4000 


1099+3 


3.5+0.2 



Table 1. Log of seven PCA observations selected for detailed 
spectral analysis. The root mean square variability is calculated 
from 3-20 keV light curves with 16-s resolution. 



numbered stars on Fig. 2, and listed in Table 1. We will refer 
to these spectra as S1-S7. The best-fitting parameters of the 
simple model Ml are shown in Table 3. We also show un- 
absorbcd luminosities of the disc (Ldisc) and Comptonized 
(icomp) components, calculated for a distance of 12.5 kpc 
and disc inclination of 66° . The disc is not significantly de- 
tected for S3-S6. The fits are generally an excellent descrip- 
tion of the data, despite their use of a phenomenological 
broad line and smeared edge to approximately model the 



reflection features. The one exception to this is S7, which 
has x^/i' = 62.2/36. Fig. 3 (panel Ml) shows the residuals 

of this fit in terms of the ratio of data/model. 

In the next step we replace the phenomenological de- 
scription of the reflection features with the reflection code 
of Zycki, Done & Smith (1998). This calculates the self- 
consistent line emission and reflected continuum from an 
ionized, relativistically smeared disc illuminated by a Comp- 
tonized continuum. We fix the iron abundance at 3 times so- 
lar (Lee et al. 2002). The model is parameterized by ampli- 
tude of reflection (0 < Q,/2-k < 2), reflector ionization, ^, and 
inner disc radius (6 < Rin/Rg < 1000, where Rg = GM/c^). 
We call this model M2 (see Table 2) . Fitting results in Ta- 
bic 4 show that the fit is generally adequate except for S7 
again, but the is always worse than for the simple model 
Ml fits. The physical description of the line and reflection 
features allows less freedom to fit the spectrum than the 
phenomenological components, and shows that there is more 
spectral complexity in the data (predominantly in S7) than 
is described by reflection. 

Fig. 3 (panel M2) shows the residuals for the refiec- 
tion model of S7, where the main mismatch is a marked 
excess in the data at ~5.5 keV. This feature is systemati- 
cally present in the other 6 spectra also, but is much less 
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Obs. fcTdisc r feTcomp n/2-7T log{^) Ri 

n -^disc -^'comp 

(kcV) (kcV) (ergcms-i) (Rg) (10^** erg g-^) 



1 


< 0.62 


r, r,7+0 02 


> 18 


0.15 ±0.04 


Q 7+0.2 
' -0.5 


43+80 


< 1.8 


c q+0.9 
"'■■'-0.6 


31.1/37 


2 


0.91 ±0.09 


2.30 ±0.02 


> 17 


17+0" 
'-'•^'-0.04 


^■^-0.7 


43+75 

^•^-23 


3g+0.15 
'-'•'-'"-0.13 


4 O+O-" 


28.9/37 


3 


1.96 ±0.25 




4 9+0.6 


> 1.26 


< 2.6 


41+40 
^^-14 






21.0/37 


4 






6.4 ±0.7 


r, 1O+0.04 

0.1S_(, Q2 


o 4+0.3 


> 100 


< 4.5 




46.6/37 


5 






> 6.2 


0.18t°°^ 


3 4+0.3 


> 100 


< 5.6 


5 6+3-3 


34.5/37 


6 




■-'■^•-'-0.07 


q cq+O.ll 


0.50«:1? 


< 1.9 


40+230 

^"-24 


< 1.9 


10 7+^-0 


16.8/37 


7 


1 -17+0.03 
^•^'-0.06 


1 c+O.S 
^•^-0.3 


2.58«;?« 


1.40«i«, 


3 g+O.l 
•^■"-0.2 


> 350 


6 8+°-2 


n 5+0-3 


208/37 



Table 4. Best- fitting parameters for the selected spectra using a model with relativistically smeared, ionized reflection (M2 in Table 2). 
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Figure 3. Residuals (data/model) for spectrum 7 (Table 1) fit 
with different models as summarized in Table 2. M4 is the most 
complex, but is the only physically motivated model which can 
fit this spectrum. 



significant. While such features are obviously indicative of 
a strongly redshifted, broad shoulder to the iron line emis- 
sion, these effects are already included in our model spectra 
(which has iron at 3x solar abundance, so the line is ex- 
tremely strong). Instead this most probably represents the 
effects of excess absorption. Individual edges from neutral 
elements can be resolved with the Chandra grating spec- 
tra, and these clearly show that the neutral absorption has 
strong abundance anomalies, with Si and Fe being a fac- 
tor 3.5 and 2.3 overabundant, respectively (Lee et al. 2002). 
This most probably arises from material ejected during the 
formation of the black hole (Lee et al. 2002) , so it should be 
fairly far from the system, and constant with time. We use 
the VARABS code to model this, fixing the column of all the 
elements at 4.7 x lO'^^ cm"-^ except for Si and Fe which are 
fixed at 16.4 and 10.9 x 10^^ cm~^ respectively as measured 
by Chandra (Lee et al. 2002). We call this model M3 (Table 
2). 

We do not tabulate the detailed results of this inter- 
mediate model and only note that the increased absorption 



improves the fit for no increase in the number of degrees of 
freedom for all the data except S3. While /'^ ~ 34/37 for 
S3 is somewhat worse than that for the standard abundance 
model, it is an adequate fit to the data. However, S7 still has 
an unacceptable ^ jv = 93/37, as shown by its residuals in 
Fig. 3 (panel M3). 

High- and moderate-resolution spectral data indicate 
yet more complexity in the absorption, in that as well as 
the anomalous abundance neutral absorption there is clear 
evidence for ionized absorption. Both ASCA and Chandra 
detect resonance line absorption from highly ionized iron, 
and show that this is variable with time and/or spectral state 
of GRS 1915+105 (Kotani et al. 2000; Lee et al. 2002). We 
model this by narrow Gaussian lines corresponding to He- 
and H-like iron, with Ka transitions at 6.67 and 6.9 keV, 
respectively. We also include the (fixed energy) K/3 and 7 
lines at 7.80, 8.21 keV for He-like and 8.16, 8.61 keV for H- 
like, with fixed relative intensity to the Ka transition of 0.07 
and 0.06 for He-like and 0.2 and 0.07 for H-Uke (Kotani et al. 
2000). Thus there are only two additional free parameters, 
which are the intensities of the 6.67 and 6.90 keV absorption 
lines, denoted as 76.67 and /6.90, respectively. These are not 
strongly detected in Sl-6, but the fit to S7 is significantly 
improved, and now gives an acceptable 'x^ /v = 34.2/35. 
Residuals to this fit are shown in Fig. 3 (panel M4). 

We use the xstar code (Kallman et al. 1996) to check 
the assumption that the dominant opacity is from the iron 
absorption lines (as opposed to lines from other elements 
and/or photoelectric edges). We assume the same Si/Fe 
abundance anomalies as in the neutral absorption, illumi- 
nated with a r = 2 power law. This self-consistent absorp- 
tion model gives an identical to that of the resonance 
absorption lines model of Table 5. The resulting parameters 
imply these lines arise in an equivalent hydrogen column of 
^ 3 X 10^3 cm~^ with ionization parameter ^ ^ lO^'^ ergs 
cm s~^. None of the seven spectra are significantly better fit 
with the full XSTAR ionized absorption models so we use the 
absorption lines in all subsequent fits as they are easier to 
constrain. 

Our final model M4 consists of the disc emission, 
its thermal Comptonization, ionized and relativistically 
smeared refiection, cold absorption with anomalous abun- 
dance of silicon and iron and resonance absorption lines from 
highly ionized iron. The fit results to all the seven spectra 
are in Table 5. 

The effective increase in cold absorption affected mostly 
the soft part of the spectrum, decreasing the soft colour by 
~0.2, when M2 was replaced by M4. The extrapolated bolo- 
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metric disc luminosity is also heavily dependent on the as- 
sumed absorption model. With model M4 the disc luminos- 
ity is significantly larger than with M2. This is particularly 
well seen for spectrum S3, where changing the absorption 
model has switched the fit from a high disc temperature 
(w 2 keV), low Comptonized luminosity model, to a low 
disc temperature (« 0.7 keV), high Comptonized luminos- 
ity model. In S6 the best-fitting M4 model is disc-dominated 
(-^0.1 of total power in the Comptonized tail), though this 
particular fit is very poorly constrained and a more Comp- 
tonized spectrum (~0.6 luminosity in the tail) is also possi- 
ble. This shows the difficulties in dealing with such an ab- 
sorbed source, and illustrates the reason for the conservative 
limit of DG03 of restricting their sample to objects with col- 
umn of < 2 X 10^^ cm~^. 

We note that none of the models including reflection 
(M2-4) give any indication for the extreme rclativistically 
smeared components expected for an X-ray illuminated disc 
which extends down to the last stable orbit around a black 
hole. All the inner radii derived from relativistic smearing 
arc substantially larger than -^20 Rg (10 Schwarzchild radii). 
Similarly large radii are also inferred by Martocchia et al. 
(2002) from BeppoSAX data on this source. This is impor- 
tant as there is some evidence that GRS 1915-fl05 is an 
extreme Kerr black hole (Zhang et al. 1997; Cui et al. 1998; 
Sobolewska & Zycki 2003), with the disc extending down 
to 1.24 Rg. However, the line broadening does not rule out 
such a disc since the reflector is highly ionized. The inner- 
most disc may well be so highly ionized that iron in the disc 
is completely stripped, so there are no line features produced 
in the most strongly curved space-time regions though there 
can be continuum emission/reflection from the disc (see also 
Wilson & Done 2001 for a similar conclusion for the very 
high state spectra of XTE J1550-564). 

3.4 Physical model applied to all pointed spectra 

The previous section has shown that a physically motivated 
model for GRS 1915-f 105 which can adequately describe all 
the selected spectra includes anomalous abundance neutral 
absorption, ionized absorption, relativistically smeared, ion- 
ized refiection, with continuum from the disc and thermal 
Comptonization. The high absorption means that the de- 
rived soft colour and disc luminosity change substantially 
with changing the model. Nonetheless, we now have a phys- 
ical model which incorporates all the observed spectral fea- 
tures so we can at least use this to get our best estimate for 
the intrinsic colours of the system. We refit the 292 spec- 
tra from Fig. 2 with this physical model, and the resulting 
colour-colour diagram is shown in Fig. 4. Black and blue 
points on this diagram correspond to spectra which show 
little variability, so the model parameters represent the phys- 
ical components of the spectra. However, this is not the case 
for the spectra averaged over high- variability periods (shown 
in green and red in Figs. 2 and 4), where the physical model 
is merely parameterizing the time averaged spectrum to give 
a best estimate for the colours. A comparison with Fig. 2 
confirms the result from the individual spectral fits that the 
major change is that the soft colour is reduced by ~0.2 in 
all spectra, so that the overall pattern stays the same. 

Fig. 4 also shows the best-fitting unfolded individual 
spectra, together with the model components. Spectra along 



the upper diagonal branch between SI and S2 have a hard 
component which has no discernable rollover in the PCA 
bandpass (fcTcomp 10 keV), with spectral index steepen- 
ing from SI to S2. S3-S5 form a parallel diagonal track 
for changing spectral index, but here the electron temper- 
ature is clearly seen in the high-energy spectral curvature 
(fcTcomp ~ 6 keV). All these spectra have comparable disc 
and hard tail luminosity. S6 and S7 are different from S1-S5 
in that there is an increase in disc temperature. The hard 
tail is characterized by even lower temperature electrons, 
kTcomp ~ 2-3 keV. The disc temperature > 1 keV rather 
than hard tail dominates the change in soft colour. 

Only spectra with colours similar to those of S7 show a 
strong detection of the resonance absorption lines from He- 
and H-like Fe. This is consistent with previous moderate- 
and high-resolution spectra from ASCA and Chandra. The 
Chandra spectrum of Lee et al. (2002) has a continuum simi- 
lar to S2, with resonant line intensities of a few 10~^ photons 
cm~^ s~^, easily consistent with our upper limit in Table 5. 
The three ASCA spectra of Kotani et al. (2000) taken in 
1994 September, 1995 April and 1996 October, have colours 
similar to S7, midway between S6 and S7, and similar to 
S5, respectively. The lines are only significantly detected in 
the 1994 and 1995 ASCA data, where they have equivalent 
widths of 30-50 eV for the He- and H-like Ka lines. This is 
somewhat smaller than, but comparable to the 100-120 eV 
equivalent widths for these lines we derive for S7. This ma- 
terial is almost certainly part of an outflowing wind, which 
can have a comparable mass loss rate to that of the mass 
accretion rate required to power the observed emission (Lee 
ct al. 2002). That the lines are most visible in S7 can be 
explained as an ionization effect. All the other spectra are 
harder, so a higher proportion of iron can be completely ion- 
ized, giving smaller equivalent width features. This can also 
feedback onto the wind acceleration mechanism if there is 
considerable line driving, so that there is intrinsically less 
wind at higher ionization (e.g. Proga & Kallman 2002). 



4 GRS 1915-1-105 IN CONTEXT 

While our physically based models give generally good fits 
to all the individual observation ID spectra, these data are 
extracted from full RXTE pointings. These have exposure 
times of several kiloseconds, so many of them average over 
periods of strong spectral variability. Thus the parameters 
derived from the spectral fits in Fig. 4 will not show the ex- 
tremes of luminosity, and can be misleading since the model 
components are not linear. 

Thus we need to go back to the 128-s spectra in order 
to get a clearer idea of the properties of GRS 1915-1-105. 
However, it is not feasible to fit all 16489 spectra with the 
computationally intensive physical model, where a single fit 
can take tens of minutes. Fortunately, the spectral fits from 
Section 3.3 show that most of the change in derived soft 
colour and soft luminosity between the simple phenomeno- 
logical description and the more physical model arises from 
the change in cold absorption. Yet most of the increase in fit- 
ting time comes from modelling reflection. Hence we use the 
simple phenomenological model of smeared edge and broad 
line to describe reflection (as in model Ml of Table 2), with 
fixed anomalous abundances of Si and Fe to test how well 
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Table 5. Best-fitting parameters for the selected spectra (Table 1) using our best physical model -with reflection, anomalous neutral 
abundances and ionized iron absorption lines (M4 in Table 2). An asterisk indicates that the parameter is completely unconstrained. 



this can reproduce the colours and luminosities derived from 
the physical model (we call this model M5; see Table 2). A 
fit of this simple model to the seven selected spectra from 
Section 3.3 shows that it gives a very good match to the 
colours from the physical model, and reproduces the soft lu- 
minosities to within a factor of ~2. This gives some measure 
of the systematic uncertainty in deriving the total (model 
dependent) luminosity. 

We use this simplified model (M5) to fit all the 128-s 
spectra. The left panel in Fig. 5 shows the colour-colour plot 
derived from this (black, red and yellow points), compared 
to the colours derived from the 'normal' black holes of DG03 
and GD03 (cyan points). DG03 show how position on this di- 
agram relates to the different black hole spectral states (see 
e.g. the reviews by Tanaka & Lewin 1995; Esin, Mclintock & 
Narayan 1997), which wc briefly summarize here. The power 
law dominated low/hard state corresponds to the diagonal 
track starting from colours (2.0, 1.5) for an photon spec- 
tral index of F ~ 1.5, extending down to about (1.5, 1.0) 
as the spectrum softens to F 2. Intermediate or very high 
state spectra have roughly equal luminosity in a disc com- 
ponent as in a fairly steep (F ~ 2-2.5) power-law tail and 
have colours within ~0.2 of (1.3, 0.8). The disc is completely 
dominant in the ultrasoft state, so the colours are close to 
those predicted by the pure disc blackbody track. The clas- 
sic high/soft state (as seen e.g. in Cyg X-1: Gierlinski et al. 
1999) has a disc-dominated soft spectrum but at higher en- 
ergies has a noticeable fraction of luminosity in a power-law 
tail. This gives it somewhat harder colours than the ultra- 
soft states at ~(0.7, 0.8). This is the region identified by 
DG03 as a sufficient (but not necessary) condition for the 
object to be a black hole. 

The spectra from GRS 1915+105 with colours similar to 
S1-S5 overlap with the intermediate/very high state region 
on the colour-colour diagram, while S6 and S7 occupy the 
high-temperature end of the ultrasoft region. GRS 1915+105 
never goes into the classic low/hard state ^(2.0, 1.5), nor 
does it go into the 'black hole only' region of DG03 (the 
classic high/soft state). 

We categorize each 128-s spectra according to variabil- 
ity of a period they come from, marking them in black, red 
and yellow for r.ni.s. < 40, 40-60 and > 60 per cent, respec- 
tively. This makes it clear that the major instability track is 
between spectra with colours like S2 and S6 (see Fig. 2). A 
comparison with Fig. 4 shows that the higher time resolution 
of these spectra (128-s as opposed to a few ks) gives a simi- 
lar diagonal track for the limit cycle, but extending further 



towards the disc blackbody line (i.e. corresponding to higher 
temperatures and luminosities) . With even higher resolution 

data then this track extends even further (sec the solid and 
dashed line in Fig. la), showing that our 128-s spectra do 
not follow the fastest variability of GRS 1915+105. 

The upper-left end of the track corresponds to spectral 
state C (Fig. lb), rather mislcadingly also called the hard 
state. In fact, it is very similar to e.g. very high state data 
from XTE J1550-564 (Gierlinski & Done 2003a). There is 
nothing special in the spectral shape here which hints at 
the origin of the unique variability of GRS 1915+105. The 
other end of the instability track lays in the ultrasoft re- 
gion, corresponding predominantly to the higher temperar- 
ture state B (Fig. lb). These have similar colours to pure 
disc spectra, but are at a higher effective temperature than 
any of the other black holes. Thus the only hint of differ- 
ence in the colours for the origin of the dramatic instability 
in GRS 1915+105 is that the ultrasoft state appears to ex- 
tend to higher disc temperatures. 

The right panel of Fig. 5 is a colour-luminosity plot for 
GRS 1915+105 superimposed on that for all the other black 
holes (DG03, GD03). The colour-coding is the same as in the 
left panel. This is much more illuminating as to the origin of 
the unique instability behaviour of GRS 1915+105. Here it is 
plain that the instability occurs in the most luminous spec- 
tra, those which have L > Z/Edd, and only GRS 1915+105 
goes up to such high luminosities. The luminosity never 
drops below ~0.3 LEdd, explaining why GRS 1915+105 
never goes into the hard state (typically seen only at lu- 
minosities lower than a few percent of Eddington). It also 
explains why GRS 1915+105 never goes into the 'black 
hole only' area of the colour-colour diagram. Spectra in 
this region, characterized by a low-temperature disc com- 
ponent together with a weak hard tail, are seen predomi- 
nantly around ~0.1 I/Edd- The persistently high luminosity 
of GRS 1915+105 takes it above these low disc tempera- 
tures. 

Where GRS 1915+105 has a luminosity which over- 
laps with that of the 'normal' black holes it shows com- 
parable spectra and does not show the dramatic limit cycle 
variability. However, it also extends to higher luminosities, 
up to ~2 3 I/Edd- This makes it unique among the black 
holes studied by DG03 and GD03, and it is here that the 
unique instability occurs. However, supereddington luminos- 
ity seems to be only a necessary, not sufficient condition for 
the disc instability. There are several low variability spec- 
tra from GRS 1915+105 which exceed the Eddington limit. 
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Figure 4. The intrinsic colour-colour plot for all tlic observations sliown in Fig. 2 fitted with the physical model including reflection, 
anomalous neutral abundances and ionized iron absorption lines (M4, Table 2). The individual spectra corresponding to the selected 
observations listed in Table 1 are shown deconvolved with this model, where the red, blue and green lines give the disc, Comptonized 
emission and its reflection, respectively. The best-fitting parameters for each spectrum are shown in Table 5. 



Perhaps the instability requires that the disc as opposed 
to total luminosity goes above Eddington, or perhaps this 
merely shows that the disc structure takes some time to re- 
spond to an increase in central mass accretion rate (e.g. van 
der Klis 2001). 

The left panel in Fig. 5 includes the background colour 
tracks for a disc blackbody spectrum, and the high luminos- 
ity, ultrasoft spectra on the extreme end of the limit-cycle 
variability intersect this line at colours corresponding to 
fcTdisc ^ 2.2 keV. A more detailed analysis is complicated as 
many of these most extreme colour 128-s spectra still average 
over substantial variability. However, these ultrasoft spectra 
have similar colours to the stable spectrum S6. // they are 



truly similar then the results in Table 5 indicate that these 
are not simply high-temperature disc emission. Instead they 
might be much better fit by a lower temperature disc at 
< 1.6 keV, with a fairly strong ( ^ 10 per cent of bolometric 
flux), low temperature (few keV) Comptonized tail (see also 
Zdziarski et al. 2001). A similar effect is also seen in the 
high-temperature ultrasoft spectra of GRO J1655-40 (Kub- 
ota, Makishima & Ebisawa 2001; GD03). 

However, irrespective of the detailed spectral form of 
the highest luminosity spectra, we include the minimum ab- 
sorption column, so derive (on average) a lower limit to the 
flux. Additional, time variable cold absorption components 
(e.g. Sobolewska & Zycki 2003) will only increase the in- 
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Figure 5. The intrinsic colour-colour (left) and colour-luminosity (right) diagrams for GRS 1915+105 created from 128-s spectra. Black, 
red and yellow points denote periods of low, high and very high variability with r.m.s. < 40, 40-60 and > 60 per cent, respectively. On top 
of this we have superimposed the data from 'normal' black holes (cyan points) from DG03 and GD03. The blue lines on the left diagram 
show the model colours of a power law (upper line, with F increasing from 1.5 to 3.0 diagonally downwards) and a multicolour disc (lower 
line, with fcTjiat. increasing from 1 to 3 keV rightwards). The yellow points track the highest variability of GRS 1915-1-105, related to 
limit-cycle instability. It occurs between colours corresponding to the very high state (1.1, 0.8) and a hot ultrasoft state corresponding 
to an effective disc temperature of excess of 2 keV. 



ferred luminosity. There is a factor ~2 systematic uncer- 
tainty on the distance (7-12 kpc: Fender et al. 1999, though 
the lower limit could present significant problems for the de- 
rived luminosity of the companion star: Greiner et al 2001b) 
so this could reduce the bolometric luminosity estimates by 
at most a factor 4. Thus it seems that the trigger for the 
characteristic limit-cycle variability in GRS 1915-1-105 is al- 
most certainly super-Eddington accretion rates. 



5 LIMIT-CYCLE VARIABILITY 

GRS 1915-H05 is the most luminous black hole system stud- 
ied here, and is also the only one to show the characteristic 
limit-cycle variability. Plainly the results are consistent with 
the accretion flow onto a black hole becoming unstable at 
very high luminosities of > I/Edd- The same instability does 
not seem to operate in neutron stars. Z sources emit at 1-3 
^Edd, while Cir X-1 can reach up to lOI/Edd (DG03). Ap- 
proximately half of the neutron star luminosity should be 
from the boundary layer, so the accretion flow luminosity 
is only ~0.5-1.5 I/sdd in the Z sources, perhaps too low to 
trigger the instabihty (but not in Cir X-1). Alternatively, 
irradiation of the accretion flow by the boundary layer may 
act as a stabilizing mechanism (Czerny, Czerny & Grindlay 
1986, but see Mineshige & Kusenose 1993 and below). 

Following Belloni et al. (1997a, b), the instability it- 
self is generally associated with the well known radiation 
pressure instability of a standard a disc at high luminosi- 
ties (Lightman & Eardley 1974; Shakura & Sunyaev 1976). 
In its simplest form this predicts that the inner accre- 



tion flow cycles between a highly luminous, hot, advection- 
dominated disc and a much lower luminosity, cooler, gas 
pressure-dominated disc (Abramowicz et al. 1988; Honma, 
Matsumoto & Kato 1991; Szuszkiewicz & Miller 1997, 1998; 
Zampieri, TuroUa & Szuszkiewicz 2001). One key problem 
with this interpretation is that the limit-cycle variability 
should occur at well below I/Edd (see Nayakshin, Rappaport 
& Melia 2000 for other problems). Theoretical models show 
that the limit cycles become noticeable at Ldiac ^ 0.3I/Edd 
(Janiuk, Czerny & Siemiginowska 2002), yet none of the 
comparison sample of 'normal' black holes above this lumi- 
nosity show the limit cycle (Kubota et al. 2001; Kubota & 
Makishima 2003; GD03). 

GD03 summarize various ways to delay the onset of the 
radiation pressure instability in a standard disc. Generally 
these involve some additional energy loss channel e.g. dissi- 
pating some of the accretion energy in a hard X-ray corona, 
or via a wind or jet. The problem with these is that the 
'normal' black holes can show disc-dominated spectra (so 
no obvious energy loss in a corona/jet) without producing 
the instability even at Ldisc ~ 0.5Z/Edd (GD03). By contrast 
GRS 19154-105 shows strong coronal emission and has a jet 
and a wind, all of which should act to stabilize the disc, 
yet does show the instability. Irradiation of the disc by the 
strong coronal emission should also act as a stabilizing mech- 
anism, similarly to that proposed for the boundary layer in 
neutron stars (Czerny et al. 1986). Thus it seems unlikely 
that irradiation alone can stabilize the disc in the Z sources 
since it does not manage to stabilize GRS 1915-1-105. 

The 'normal' black holes show that the disc is stable up 
to at least Ldisc ~ 0.5LEdd, so its structure must be rather 
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different from that predicted by the standard a disc equa- 
tions (GD03). The fact that the instabihty in GRS 1915+105 
is triggered at or near Eddington luminosity seems to sug- 
gest that radiation pressure does play a pivotal role in the 
formation of the dramatic variability. That the a disc equa- 
tions do not work in detail is perhaps not surprising given 
the ad hoc nature of viscous heating in these models. The 
next generation of disc models will give a much better es- 
timate for the disc structure, as these can incorporate the 
MHD dynamo which is the physical origin of the viscosity 
(Balbus & Hawlcy 1991). When these models are sufBciently 
developed to solve the time dependent 3D, coupled radia- 
tion and magneto-hydrodynamic equations over a substan- 
tial range of radii (Turner, Stone & Sano 2002; Turner et 
al. 2003) then the question of inner disc stability for highly 
luminous systems can be revisited. 



6 LONG-TERM VARIABILITY 

The unique short-term limit-cycle 

variability of GRS 1915-1-105 is plausibly because it reaches 
higher L/L^dd than any other known Galactic black hole. 
But then the question becomes why GRS 1915-f 105 alone 
reaches such high luminosities? The answer to this lies in 
the evolutionary state of GRS 1915-1-105. The system has 
the longest known period of any LMXB of 33.5 days, imply- 
ing a very wide orbit (Greiner et al. 2001a). The secondary 
is a low-mass giant star (Greiner et al. 2001b), so the ac- 
cretion takes place via Roche lobe overflow (Eikenberry & 
Bandyopadhyay 2000). 

There is now considerable agreement that transient out- 
bursts in the LMXBs are caused by the classic disc instabil- 
ity mechanism which operates when hydrogen goes from be- 
ing mainly neutral to mainly ionized. This is a very different 
physical mechanism for a disc instability than the radiation 
pressure instability discussed above, and happens at much 
lower luminosity, but it gives rise to another S-shaped ther- 
mal equilibrium curve (Smak 1982). In discs around white 
dwarfs it can trigger limit cycles where there are disc out- 
bursts on timescales of days-weeks, followed by a period of 
quiescence lasting weeks-months (see e.g. the review by Os- 
aki 1996). In LMXBs the disc structure is strongly modified 
by irradiation once the outburst starts, keeping hydrogen 
ionized for much longer so that the outburst timescales are 
dramatically extended (van Paradijs 1996; King 1997; King 
& Ritter 1998; Dubus et al. 2001). 

The peak luminosity at the start of the LMXB outburst 
is determined by the disc size, Rt- The maximum mass the 
disc can hold before going into outburst is oc Rt, while the 
inflow timescale is oc Rt- Thus the peak luminosity oc Rt 
so the bigger the disc, the bigger the peak luminosity (King 
& Ritter 1998). But the disc size in LMXB is fixed by tidal 
torques to ~1.4x_Rcirc, where -Rcirc is the circularization ra- 
dius of the binary (Shahbaz, Charles & King 1998), and the 
binary system parameters are determined by the require- 
ment that the sccondarjr star can overflow its Roche lobe. 
For main sequence secondaries the orbit must be fairly small, 
so the disc is small. GRS 1915-1-105 has a giant secondary, 
so has a huge disc, with Rt ~ 10^^ cm for a mass ratio 
of 0.1 (e.g. Frank, King & Raine 1992). The quiescent disc 
mass is ~10^* g (Shahbaz, Charles & King 1998), which can 



maintain an Eddington accretion rate for ~10 years, similar 
to observed outburst timescale. 

Thus both the long and short-term unique variability 
behaviour of GRS 1915-1-105 can be explained as a function 
of its evolutionary state, where the wide orbit allows a huge 
disc to form. 



7 CONCLUSIONS 

GRS 1915-1-105 is not in a separate class from 'normal' black 
holes. When it is at the same L/L-^dd then its spectra and 
(more importantly) time variability behaviour are similar to 
that seen in the 'normal' black holes. Its unique limit cycle 
variability only appears when it radiates at uniquely high 
(super Eddington) luminosities. 

GRS 1915-1-105 boldly goes to luminosities where no 
black hole has gone before, and we eagerly await the next 
generation of disc models which will more clearly identify 
the origin of the instability. 
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